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Abstract: Considering the aging effects of existing Inner Drift Chamber (IDC) of BESIII, a GEM based inner 
tracker is proposed to be designed and constructed as an upgrade candidate for IDC. This paper introduces a full 
simulation package of CGEM-IT with a simplified digitization model, describes the development of the softwares for 
cluster reconstruction and track htting algorithm based on Kalman hlter method for CGEM-IT. Preliminary results 
from the reconstruction algorithms are obtained using a Monte Carlo sample of single muon events in CGEM-IT. 
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1 Introduction 

The Beijing Spectrometer HI (BESIII) [1] is a de¬ 
tector of high precision that operates at the Beijing 
Electron-Positron Collider II (BEPCII) [2], a high lu¬ 
minosity, multi-bunch e+ e“ collider running at the r- 
charm energy region. The Main Drift Chamber (MDC) 
of BESIII is a key subdetector for charged track recon¬ 
struction with good spatial and momentum resolutions. 
MDC comprises two parts, Inner Drift Chamber (IDC) 
and Outer Drift Chamber (ODC). However, the existing 
IDC is suffering from aging problems after five years of 
operation, which results in an obvious decrease in its gas 
gain and spatial resolution [3]. 

Gas Electron Multiplier (GEM) detector [4] is char¬ 
acterized by a significant radiation tolerance, high rate 
capability, an outstanding spatial resolution and a lim¬ 
ited radiation length so that it can serve as an ideal 
inner tracker. Then a Cylindrical-GEM Inner Tracker 
(CGEM-IT) is proposed to be an option for the upgrade 
of the IDC for BESIII. The proposed CGEM-IT con¬ 
sists of three independent layers of CGEM detectors with 
outer radii 87.5mm, 132.5mm and 175mm and lengths 
532mm, 690mm and 847mm in axial direction. Each 
CGEM detector is composed by a cathode, three GEM 


foils and a two-dimensional readout anode in a concen¬ 
tric arrangement, as shown in Figure High voltages 
are applied on each GEM foil full of micro bi-conical 
holes, providing strong electric fields in the holes to am¬ 
plify the ionized electrons from a track to be measured. 
Such an inner tracker combined with ODC is expected 
to offer excellent spatial resolutions: about 100 /im in 
the r0 view and about 200 jam in the z direction, where 
r is the radius, 0 is the azimuth angle and z is the coor¬ 
dinate along the central axis of a CGEM detector. The 
development of simulation and reconstruction softwares 
for CGEM-IT is essential to study its performance in 
track measurement and used for design optimization of 
the detector. 

In this paper we will introduce the full simulation 
package for CGEM-IT, describe the implementation of 
the cluster reconstruction and track fitting algorithms 
for this detector in the BESIII Offline Software System 
(BOSS) [5] framework, and present the preliminary re¬ 
sults on the resolutions of CGEM-IT and the perfor¬ 
mance of reconstructed tracks by using a Monte Carlo 
sample of single-muon events. 
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Read out 



Fig. 1. The structure of one CGEM detector that 
consists of one cathode, three GEM foils and one 
anode. When a track (red solid curve) traverses 
one CGEM, ionization occurs in the drift area (the 
gap between the cathode and the hrst GEM foil). 
The electrons are multiplied in the holes of GEM 
foils and drift towards the readout anode. The 
position in anode plane is obtained by digitiza¬ 
tion with (black dashed lines) and without (blue 
solid lines) taking into account the Lorentz angle 
Oi and electron diffusion. 


pie scattering and so on. The information of each step 
in drift area is recorded as Cgem Hit. These Cgem hits 
are summed up to get the position and total energy loss 
of the track segment in the drift area as Cgem Truth. 
Then an algorithm called Cgem Digitizer converts the 
Cgem Truth to the electronic readout signal named as 
Cgem Digi. 
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2 Simulation of CGEM-IT 


Eig. 2. Elow diagram for the simulation of CGEM-IT 


Reliable Monte Carlo (MC) simulation plays an im¬ 
portant role in describing the geometry and material of 
a detector and its response in track measurements. A 
full simulation package based on GEANT4 [6] has been 
developed for CGEM-IT and the BOSS framework has 
been modified correspondingly. 

2.1 Simulation procedure 

Figure 1^ shows the flow chart of the simulation proce¬ 
dure. The geometries and materials of CGEM-IT includ¬ 
ing the micro holes in GEM foils are fully and exactly 
constructed according to the preliminary design of the 
CGEM-IT. Each GEM foil is a copper clad Kapton foil 
as thin as 54 /im with holes evenly distributed on it. The 
holes are of the same size with an external diameter of 70 
/im. The space between the cathode and the first GEM 
foil is defined as the drift area, i.e. sensitive volume. The 
simulation is processed event by event. The final states 
(primary tracks) and their initial momenta of each event 
are produced by an event generator. 

To simulate the process that a charged primary track 
traverses the CGEM-IT, the track in drift area is prop¬ 
agated in discrete steps along the trajectories, during 
which all possible interactions between track and CGEM- 
IT are considered including ionization. Coulomb multi- 


2.2 Digitization 

The electrons ionized by a charged particle crossing 
the drift region will be multiplied in the micro holes of 
the GEM foil and drift towards the anode where signals 
are induced. Three GEM foils are assembled to multiply 
the number of the electrons by a factor up to 10® [3]. 

Both the magnetic field and diffusion of the electrons 
have impact on the drift properties and charge distribu¬ 
tion in the readout anode [7] . 

The presence of an axial magnetic field makes the 
drift line deviate from the direction of the electric field 
and the deviated angle is defined as Lorentz angle, as 
shown by the dashed lines in Figure The effect of 
diffusion will enlarge the cluster of multiplied electrons. 
In this work, a simplified digitization model is used, in 
which the Lorentz angle a and electron diffusion are not 
considered. As shown by the solid line in Figure the 
entrance and exit points in the drift area are projected 
to the anode plane as the extent of the signal in the read¬ 
out anode. The readout anode of each layer is segmented 
with 650 /im pitch XV patterned strips with a stereo an¬ 
gle that changes depending on the layer geometry. The 
fired strips are determined by the size of the cluster of 
multiplied electrons. The charge deposited in the drift 
area is multiplied by an expected gain and shared by the 
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fired strips. The identifier of each fired strip as well as 
its collected charge is recorded as one Cgem digi just like 
real experimental data. 

There are two different strip readout methods, the 
binary one and the analog one [3], lead to different digi¬ 
tization models. The binary readout method sets a fixed 
threshold on the total charge of the fired strips and the 
reconstructed position is the geometrical center of these 
strips. In the analog readout method, 

the distributions of charge deposited along two read¬ 
out coordinated will be recorded. The collected charge of 
each strip is obtained by the integration of the distribu¬ 
tions so that one can set a threshold on the single strips 
and reconstruct the charge centroid of fired strips as the 
position of the cluster, which results in a better spatial 
resolution compared with the binary readout method. 
The precise shape of the charge distributions will be de¬ 
termined by the beam test on the CGEM planar proto¬ 
type afterwards. So in this paper, the digitization model 
is based on the binary readout method, where total de¬ 
posited charge are shared equally by the fired strips. 


2.3 Modification of framework 


As CGEM-IT is a new detector added to BESIII, the 
BOSS framework doesn’t provide necessary interfaces for 
the softwares related to CGEM-IT. Therefore we update 
the framework by adding the components for data pro¬ 
cessing with CGEM-IT. The information of Cgem Digi 
and Cgem Turth is stored in Transient Data Store (TDS) 
and can be accessed by the reconstruction algorithms via 
the event data service. 


3 Cluster reconstruction for CGEM-IT 


When the multiplied electrons reach the readout an¬ 
ode, a group of X and V strips will be fired. One cluster 
is defined in the anode plane and consists of continu¬ 
ous fired strips in X direction and V direction that have 
intersections. Eigurej^ exhibits one example of the com¬ 
position of one cluster. Next we will discuss the cluster 
reconstruction and the obtained spatial resolution. 



Fig. 3. Schematic of one cluster induced by a 
charged track, composed of 3 X-strips and 3 V- 
strips on an unrolled anode plane, where the 
stereo angle 0 is the angle between X-strip and 
V-strip, the X-axis is perpendicular to X-strips, 
the V-axis is perpendicular to V-strips and the 
C-axis is parallel to X-strips. 

3.1 Cluster reconstruction algorithm 

Eirstly, the information of fired strips is retrieved 
from the raw data, including the strip ID, strip type 
(X or V) and collected charge. Then sets of continuous 
fired strips in X and V directions are found and defined 
as X-clusters and V-clusters respectively. Next the pos¬ 
sible intersections between X-clusters and V-clusters are 
searched for and saved as XV-clusters. The position in X 
(V) of each cluster is acquired by calculating the charge 
weighted geometrical center of the fired strips that com¬ 
pose this cluster: 



(1) 


EfEm 

(2) 

EfEQ. 


where i is the index of strips, Xi (V) is the middle po¬ 
sition of strip i, Qi is the charge collected by each strip 
and Nx {Ny) is the number of fired X-strips (V-strips). 
The corresponding position in the absolute cylindrical 
coordinates (r, 0, z) is given by 

r = R 

< z = {V — XcosO)/sinO (3) 

. <t^ = X/R 

where R is the radius of the anode in which the clus¬ 
ter is found, 0 is the stereo angle between X-strips and 
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V-strips. To obtain the middle positions where charged 
particles pass the CGEM in the drift area, a correction 
to the reconstructed positions of clusters is demanded. 
According to the simplified digitization model, the posi¬ 
tion of a track in drift area has the same azimuth angle 
(j) and the same z as the cluster, but with a different ra¬ 
dius. The corresponding output data model is defined 
by a class RecCgemCluster, including the following in¬ 
formation: layer ID, GEM ID, cluster type (X, V or XV), 
collected charge and reconstructed position. 

3.2 Spatial resolution for CGEM-IT clusters 

The spatial resolution for GGEM-IT clusters is an es¬ 
sential performance indicator of the cluster reconstruc¬ 
tion algorithm. The residuals in X and V directions are 
defined as: 

^rec ^truth ( 4 ) 

Sv ^rec ^truth ( 5 ) 

where Xrec and Kec are the reconstructed positions in 
the drift layer in X-direction and V-direction while Xtmth 
and Vfruth are the true positions. Eor a Monte Garlo 
sample of with a momentum of 1.0 GeV/c, the dis¬ 
tribution of the residuals in X and V are shown in Eigure 
m Note that the distribution of residuals in X is not a 
pure Gaussian, but rather a superposition of rectangu¬ 
lar and Gaussian. The Gaussian part (filled region in 
EigQa)) comes from clusters with multiple strips while 
the rectangular part (dashed line in EigJ^a)) comes from 
clusters with single strip. The RMS (Root Mean Square) 
of this curve is 217.1 /im and defined as the resolution in 
X-direction. A fit of the distribution of residuals in V to 
a Gaussian function gives a resolution of 129 /im. 




Fig. 4. The position residuals in X (a) and V (b) 
of the cluster for a Monte Carlo sample of sin¬ 
gle 1.0 GeV/c /i“ tracks. The contribution from 
clusters with multiple strips (hlled region) and the 
part from clusters with single strip (dashed line) 
are shown in (a). The solid line in (b) is the Gaus¬ 
sian ht curve for the residuals distribution in V- 
direction. 

4 Track fitting by Kalman filter with 
CGEM-IT clusters 

In a uniform magnetic field a charged particle moves 
along a helical trajectory. So a track can be described by 
a helix model with 5 parameters (dp, 00,/^,^^, tan A) [8], 
where dp and are the signed distance of closest point 
to the pivot hi x — y plane and z direction respectively, 
00 is the azimuth angle of the pivot to the helix center, 
n is 1/pt with the same sign as the charge of the track 
and tan A is the slope of the track. 

The track finding algorithm mis] and track fitting al¬ 
gorithm by Kalman filter uni in service using MDG hits 
can be easily adopted to accomplish the same task using 
ODG hits only. To combine the cluster measurements by 
CGEM-IT, the adopted track fitting by Kalman filter is 
extended to be able to extrapolate the track found with 
ODG hits into CGEM-IT and fit the track with these 
additional clusters in a right sequence. 

In this section, we will discuss the extension of 
Kalman filter based track fitting to use CGEM-IT clus¬ 
ters and present the obtained resolutions for the track 
parameters after the track fitting. 

4.1 Implementation of the track fitting by 
Kalman filter for CGEM-IT 

The description of geometries and materials of 
CGEM-IT is implemented in the initialization of the 
track fitting algorithm and identical with that from sim¬ 
ulation package. 

The reconstructed clusters are retrieved as input. 
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The track parameters and their covariance matrix in 
the ODC are propagated into CGEM-IT by many small 
steps, during which the effects from inhomogeneous mag¬ 
netic field, energy loss and multiple scattering in the ma¬ 
terial are included. When the track is extrapolated into 
a drift area of CGEM-IT, the measured position from 
the CGEM-IT cluster and the predicted track parame¬ 
ters are combined to update the track parameters. This 
process is repeated for the subsequent regions, with each 
iteration yielding more accurate track parameters and 
their covariance matrix, and stops at the original point. 
The obtained track parameters are taken as the recon¬ 
struction results for this track. 

4.2 Vertex resolution 

To obtain the vertex resolution of the reconstructed 
tracks, the residuals of dp and are calculated which are 
the differences between the fited values and the truth at 
the generation level. Eigure shows the distribution of 
the residuals for the Monte Carlo sample single /i“ with 
1.0 GeV/c, and a Gaussian fit yields a resolution of 192 
/im in dp and 0.33mm in d^. 



dp -do ,Tcm) 

^rec '^truth' ' 



Fig. 5. The dp (a) and dz (b) residuals for single 
1.0 GeV/c tracks with CGEM-IT. The curves 
are the results of a Gaussian ht. 


The vertex resolution in dp and d^ for single ijl~ tracks 
with different momenta for CGEM-IT and IDG is dis¬ 
played in Eigurej^ These results reflect that a tracking 
system using the CGEM-IT has a comparable dp resolu¬ 
tion and a significantly improved dz resolution relative 
to the case using IDG. 
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Fig. 6. The dp (a) and dz (b) resolutions of single 
li~ tracks with different momenta for CGEM-IT 
and IDC. 
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4.3 Momentum resolution 

The momentum resolution is represented by the dis¬ 
tribution of the momentum residual, which is defined by: 

Prec Ptruth (O) 

where Prec is the reconstructed momentum at interaction 
point and truth is the momentum at the generation level. 
The Sp distribution curve for 1.0 GeV/c muons is shown 
in Eigure as an instance. This curve is fitted by a 
Gaussian function and the obtained standard deviation 
5.5 MeV/c is regarded as the resolution of momentum. 
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5p(GeV/c) 


Fig. 7. The distribution of momentum residual 
for 1.0 GeV/c ii~ tracks is htted by a Gaussian 
function (red solid line). 

Figure summarizes momentum resolutions for sin¬ 
gle muon samples simulated with different momenta 
for CGEM-IT and IDC. These figures illustrate that 
their momentum resolutions are at a comparable level. 
CGEM-IT has a larger material budget and offers only 
three layers of measurements compared with IDC. How¬ 
ever, the good spatial resolution in two dimensions of 
CGEM-IT compensates the disadvantages so that a com¬ 
parable momentum resolution can be achieved. 
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Fig. 8. The momentum resolutions of the re¬ 
constructed muon at different momenta using full 
simulation of single events in CGEM-IT and 
IDC. 

5 Summary and outlook 

In order to solve the aging problems of IDC in BE- 
SIII, a CGEM-IT is suggested to be used to upgrade the 
IDC. In this paper, we present the implementation of the 
softwares for CGEM-IT. The general flow for CGEM-IT 
including simulation, cluster reconstruction and track fit¬ 
ting is realized. The results for a set of simulated /i“ 
events indicate that the vertex resolution in z direction 
is significantly improved by CGEM-IT and meanwhile 
the momentum resolution and vertex resolution in r di¬ 
rection stay at the same level with respect to IDC. The 
preliminary results demonstrate that the CGEM-IT and 
the implemented track reconstruction softwares can basi¬ 
cally fulfill the track reconstruction task as a replacement 
of IDC. 

The implementation of a more detailed digitization 
model based on the analog readout method will improve 
the spatial resolutions of CGEM-IT. Eurther studies 
should also be performed to reach a better understand¬ 
ing of the detector, including developing the track seg¬ 
ment finding aigorithms for CGEM-IT and joint tracking 
aigorithms for CGEM-IT and ODC, optimizing the re¬ 
construction softwares under the environment with noise 
and background. 
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